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Abstract 

After a few remarks about the problem of extracting transport coefficients from lattice QCD 
calculations, 1 report on recent developments in applying stochastic quantization and complex 
Langevin dynamics to field theories with a complex action due to a nonzero chemical potential. 
First results demonstrate that the sign problem poses no obstacle for this approach, even in the 
thermodynamic limit. I conclude with a comparison of two simple one-link models, describing a 
euclidean system at finite chemical potential and a Minkowski system in real time. 
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1. Introduction 

The euclidean lattice formulation of QCD provides a nonperturbative formulation of 
the theory of strong interactions, suitable for numerical simulations. This allows a de- 
tailed study of thermodynamic quantities, such as pressure, entropy and susceptibilities, 
and of the finite-temperature crossover between the confined and the deconfincd phase 
(see e.g. Ref. [1]). However, for some problems standard methods are not adequate. In 
this contribution I discuss two of these problems. In the first case, the questions asked 
are not easily answered from knowledge of euclidean-time correlation functions alone, 
but instead require access to particular real-time correlators such as spectral functions. 
This is relevant for transport properties of the quark-gluon plasma and for in-medium 
modifications of hadrons. In the second case the weight in the euclidean path integral 
is not real, leading to the infamous sign problem. This prohibits the use of importance 
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sampling, the framework underlying essentially all methods used to generate ensembles 
in lattice QCD simulations. This is relevant for QCD at nonzero baryon density. 



2. Transport and spectral functions 



Conserved currents, such as the energy-momentum tensor T^'' and the electromagnetic 
current j^, result in hydrodynamic behaviour on long length and timescales. Hydrody- 
namic structure is best visible in current-current correlation functions in real time, such 
as spectral functions. A nonperturbative determination of spectral functions p using eu- 
clidean correlators Ge obtained with lattice simulations, via the relation 

oo 

n I ^ [ '^'^ ( \ cosh[c^(T - l/2r)] 

J ZTT smh(w/2j) 



is quite involved. One widely used approach is the Maximal Entropy Method (MEM) [2]. 
However, in the hydrodynamical regime (w, |p| <C T) the problem is particularly difficult 
due the insensitivity of the euclidean correlator to details of the spectral function at small 
ui [3] . Moreover, the most commonly used MEM algorithm is inherently unstable at small 
energies and needs to be modified in order to be applicable [J] . There exist recent lattice 
calculations of the electrical conductivity [J] , the shear viscosity [5] and the bulk viscosity 
[5]. Transport and hydrodynamics from the lattice have been discussed in detail in two 
plenary talks at the Lattice conferences in the past two years |7l8j . Therefore I restrict 
myself to some remarks and refer to those contributions for more details. 

What euclidean correlators should be considered? Kubo relations relate transport co- 
efficients to the slope of spectral functions as the energy co is taken to zero. For example, 
the shear viscosity rj, bulk viscosity C, and conductivity a are determined by 

ri = lim C = - lim —, cr = lim ^ ^ \ (2) 

in terms of the spectral functions 
p'^''-'^-' {u) = / d4:,e-t([r^-(i,x),T«'^(0)]), = / d^x e^"*([/'(t, x), (0)]). (3) 



In analytical calculations it is sometimes useful to add an exactly conserved charge, 
such as 8 = j (fixT^^ or Q = j (Pxj^, to the operators in Eq. ([3]). For instance, for 
the bulk viscosity one may consider correlators of instead of T" [6|9j . While this is 
harmless in analytical calculations, it makes a difference when reconstructing spectral 
functions from the lattice. The reason is that lattice correlators and spectral functions 
are affected in a particular manner: if a conserved charge is added, lattice correlators will 
receive a r-independent contribution of the form TS, where S is the corresponding charge 
susceptibility, while spectral functions will receive a contribution 'E2'kljj5{lj). This singular 
contribution at the origin violates the smoothness condition assumed by most methods 
used to obtain spectral functions from euclidean correlators and should therefore be 
avoided. This issue was discussed already some time ago in the context of meson spectral 
functions and the current = ip^f^ijj [lOlllj . For massless noninteracting fermions the 
delta function contribution at the origin is equal for /3™(a;) and and cancels in 
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the combination However, interactions wiU smear out this delta function in the 

case of p" (resulting in a finite conductivity), whereas p^*^ is unchanged due to charge 
conservation (only the value of the susceptibility is affected). The cancelation therefore 
no longer holds and p^{u>) is not smooth as w ^ 0. 

As a second remark, let me note again that euclidean correlators are insensitive to 
details of spectral functions when the transport contribution is narrow and most of the 
weight is concentrated at w <C T. It is not so difficult to find cases where access to 
transport coefhcients from euclidean correlators is virtually impossible. This has been 
discussed for weakly coupled theories 0, heavy quark diffusion yL2j and, more recently, 
bulk viscosity in a weakly coupled theory or near a second order phase transition [3] . (For 
a lattice study of energy-momentum tensor correlators near the second order transition 
in SU(2) gauge theory, see Ref. [13].) From the viewpoint of the lattice, it would be 
more interesting to find instead examples where the transport peak is not narrow such 
that transport coefficients might potentially be accessible. So far this has been addressed 
mostly in theories with a gravity dual in the strong coupling limit |14ll5ll6] . Since these 
results provide an important motivation, more examples would be welcome. 



3. Finite chemical potential and the sign problem 



I now turn to euclidean field theories with a complex action due to the presence of a 
chemical potential. In this case the weight in the path integral is not real and importance 
sampling cannot be used. This is commonly referred to as the sign problem. If the weight 
is written as = |e~'^|e*''', we may consider the partition functions of the full and the 
phase-quenched (pq) theories, defined as 

^fuu = j Dct>e-^, Zpq- j D(t>\e~^\. (4) 

The phase-quenched theory has a real positive weight and is accessible with standard 
techniques. The sign problem is best demonstrated by the expression for the expectation 
value of the phase factor in the phase-quenched theory, 

(e^-)pc = 1^ = e-"^/, (5) 

where VI is the (four)volume and A/ is the difference between the free energy densities 
in the full and the phase-quenched theories. Eq. ^ shows that the average phase factor 
goes to zero in the thermodynamic limit, resulting in an overlap problem between the 
full and the phase-quenched theories. In other words, the important configurations in the 
full theory differ in an essential way from those in the phase-quenched case. The question 
is how to find those relevant configurations in a numerical simulation. 

Recently, this problem was reconsidered |17|18j using stochastic quantization and com- 
plex Langevin dynamics ^19,20.21^. This approach has a long history and was studied 
intensely in the 80's [22]. It was revived in the context of nonequilibrium quantum field 
dynamics in Minkowski spacetime in Refs. [23124125] and subsequently applied to eu- 
clidean theories at nonzero chemical potential. We have considered SU(3) gauge theory 
with heavy quarks [17 and the scalar 0(2) model 18], both at nonzero chemical poten- 
tial. 
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Fig. 1. Real part of the density in the full (left) and the phase-quenched (right) relativistic Bose gas as 
a function of chemical potential on lattices with size Af-*, with AT = 4, 6, 8, 10 (m = A = 1). 

The essence of this approach is that the weight is obtained as the equihbrium distri- 
bution of a stochastic process, described by a Langevin equation. For a scalar field the 
Langevin equation reads 

cl>(a:,e + e)^^ix,d)-e^^^+V~ev{x,e), (6) 

while in SU(iV) gauge theories, the links are updated according to 

U{e + e) = R{9)U{e), R = exp[iXa{eKa + Ve7]a)], Ka = -DaS[U]. (7) 

Here 9 = ne is the Langevin time, e is the Langevin timestep, the noise is Gaussian, 

(rK0))=O, {7^i9)r0))=26r.n', (8) 

Xa are the Gell-mann matrices and all other indices are suppressed. When the action is 
complex, the dynamics and field variables arc complexified. Explicitly, a real scalar field 
is written as </> ^ (p^ + icj)^, while an SU(A^) matrix is now an element of SL(iV, C). More 
details can be found in the references listed above. 

I now discuss the results obtained so far. I start with the most recent results in the 
relativistic Bose gas, i.e. a complex scalar field with 0(2) symmetry and a quartic in- 
teraction term, at nonzero chemical potential [18J. The action satisfies S*{^) — S{—fi). 
Just as QCD, this theory has what is known as a Silver Blaze problem [26] . At zero tem- 
perature the theory is in vacuum as long as the chemical potential fi is below the critical 
value fic (when interactions are ignored, fic = rn, the minimal energy of excitations). 
When fi > ^c, the theory enters a Bose condensed phase and the density is nonzero. 
The strict /i- independence of bulk physical quantities as long as /i < /ic, even though 
fi is present microscopically, is the Silver Blaze problem. We have analysed the full and 
the phase-quenched theory using stochastic quantization. The results for the density are 
shown in Fig. [TJ On the left, we show the density in the full theory, as a function of the 
chemical potential, for lattices of size N"^, with = 4, 6, 8, 10. The parameters in the 
scalar potential are m = A = 1, in lattice units. The inset shows a blowup of the tran- 
sition region. We observe the transition between the groundstate with zero density and 
the Bose condensed phase at /i = /ic ^1-15. A real phase transition can only occur in the 
thermodynamic limit and we find that the transition becomes sharper as the volume is 
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Fig. 2. Real part of the square of the field modulus in the full theory (left) and the average phase factor 
in the phase-quenched theory (right). 

increased, as expected. On the right, the density is shown in the phase-quenched theory. 
Here we observe a nonzero density for aU values of the chemical potential. This is similar 
to what is expected in phase-quenched QCD (when 1X1^^12 < /i < rripf/S). The cancelation 
of the /i-dependence of the density in the full theory is due to the phase factor e^'^ and, 
therefore, the result of the sign problem. We conclude that the complexity of the action 
and the Langevin dynamics deliver precisely what is expected. 

The /i-independence also hold for other observables. In Fig. [2] (left) we show the square 
of the field modulus |0p as obtained from complex Langevin dynamics. The value at 
zero chemical potential is nonzero due to quantum fluctuations. As fi is increased, the 
observable remains constant until the critical value is reached. This demonstrates the 
Silver Blaze feature in an impressive manner. On the right the average phase factor in the 
phase-quenched theory is shown. As indicated above, this observable gives an indication 
of the severeness of the sign problem: the phase factor goes to zero at larger chemical 
potential on all lattices and vanishes in the thermodynamic limit for all nonzero values 
of the chemical potential. This is precisely how the average phase factor is expected to 
behave [57]. Nevertheless, this approach can handle the sign problem without difficulties 
and there is no problem in taking the thermodynamic limit. 

Concerning QCD we have considered SU(3) gauge theory coupled to three flavours of 
quarks in the heavy quark limit, starting from Wilson fermions in the hopping expansion 
[T7|. In this limit, the complex fermion determinant is approximated as 

det Af « J]^ dct (l + /le^/^-Px) ^ det (l + /le^^/'^-Px ^) ^ , (9) 

X 

where h = (2k)^^ and Vx are the (conjugate) Polyakov loops. We emphasize that 
the gauge action is preserved completely and that the determinant satisfies the basic 
property [det M(^)]* = detM(-^). 

First results in this theory are obtained on a lattice of size 4^, with fixed /? = 5.6 and 
K = 0.12. In Fig. Owe show the (conjugate) Polyakov loops (left) and the density (right) 
as a function of chemical potential. The results indicate a transition from a confined 
low-density phase to a deconfined high-density phase. Again we address the severeness 
of the sign problem by studying the phase of the determinant. In this case we present 
the phase factor e^"^ — det A/(/i)/ det Af(— /i) in the full theory, and not in the phase- 
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Fig. 3. Real part of the (conjugate) Polyakov loops (P) and (P ^) (left) and the density (n) (right) as 
a function of fi for QCD in the heavy quark limit, with (3 = 5.6, k = 0.12, and Nf = 3 on a 4'* lattice. 
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Re g^"^ Lano:evin iteration 

Fig. 4. Left: Scatter plot of e^^*^ = det M(^)/ det Af (-/i) during the Lan gevin evolution for various 
values of fi. Right: Deviation from SU(3), Trt/Jt/4/3 during the Langevin evolution, for /i = 0.5 and 0.9. 

quenched theory. A scatter plot of the real and imaginary parts of this observable is 
shown in Fig. 2] (left). When the chemical potential is zero, the phase factor is equal to 1. 
At larger chemical potential, we find that phase fluctuations increase dramatically. Yet, 
observables such as the Polyakov loop and the density are under control, with reasonable 
errors. During the complex Langevin evolution, link variables are no longer in SU(3), 
instead they take values in SL(3, C). This can analysed by computing Tr [/|[/4/3, which 
equals 1 in SU(3) and is > 1 in SL(3, C). We show this quantity as a function of Langevin 
time in Fig. [4] (right). We observe that after the initial thermalization stage it is clearly 
distinct from 1, but remains bounded. This indicates that the sequence of configurations 
generated with complex Langevin dynamics differs in an essential way from configurations 
constrained to be in SU(3). For more discussion, see Refs. [17128) . 

4. Chemical potential versus Minkov^rski dynamics 

At zero chemical potential the euclidean action is real. In this case one can use real 
Langevin evolution and apply standard proofs to demonstrate that the method will 
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Fig. 5. Smallest nonzero eigenvalue of the complex Fokker-Planck operator in the one-link U(l) model 
as a function of /3 for various values of at k = 1/2 in the euclidean case (left) and for various values of 
the reweighting parameter p in the Minkowski case (right). 

converge to the correct distribution . Zero chemical potential provides therefore a 
useful reference point, see e.g. Fig. [2] (left). In the case of Minkowski (real-time) dynamics, 
where the weight in the path integral is e*"^, such a reference point is absent. Heuristically, 
one may therefore expect real time to be more difficult than nonzero chemical potential. 
To quantify this, we compared two simple one-link U(l) models [17j . building on the 
work in Ref. [25j . The partition functions are 

7r 

/dr 
e^"°''"[l-f KC0s(x-i//)], (10) 

— TT 

at nonzero /i and 

7r 

/dr 
___ cos a;-t-ipa: (H) 
27r 

in real time. The term px, with p integer, is a reweighting term, used to stabilize the 
Langevin dynamics |25| . In the euclidean case the Fokker-Planck equation for the complex 
distribution P(x, 9) reads 

The zero eigenvalue, corresponding to the stationary distribution e~'^, always exists. The 
solution of the Fokker-Planck equation can therefore be written as 



e 



-s 



P{x,e) = — + l^e-^"P^{x). (13) 

For the Minkowski case, —S is replaced with iS in the equations above. Convergence of 
the complex Fokker-Planck equation is determined by the eigenvalues A, which we studied 
numerically [17] . The smallest nonzero eigenvalue is shown in Fig. [5] as a function of (3 
in the euclidean case for various values of ^ (left) and the Minkowski case for various 
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values of p (right). In the first case, one sees that the eigenvalues are strictly positive, 
for all values of [3 and /i, indicating convergence of the complex Fokker-Planck equation. 
In the second case, eigenvalues turn negative for larger values of /3; the value of (3 where 
this occurs depends on the parameter p. While these results do not prove or disprove 
convergence for the complexified Langevin dynamics (for this the real Fokker-Planck 
equation for the real distribution p{x, y, 9) has to be studied), they are an indication of 
the difference between the two theories when applying complex Langevin dynamics. 
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